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The crystal of benzoic acid is comprised of tautomeric centrosymmetric dimers linked 
p • through bistable hydrogen bonds. Statistical disorder of the bonding protons is ex- 

cluded by neutron diffraction from 6 K to 293 K. In addition to diffraction data, vibra- 
r-| ■ tional spectra and relaxation rates measured with solid-state-NMR and quasi-elastic 

o ■ 

^ . neutron scattering are consistent with wave-like, rather than particle-like protons. 

O 

• ^ I We present a macroscopic-scale quantum theory for the bonding protons represented 

>>■ 

r~| ' by a periodic lattice of fermions. The adiabatic separation, the exclusion principle, 

and the antisymmetry postulate yield a static lattice-state immune to decoherence. 
According to the theory of quantum measurements, vibrational spectroscopy and re- 
laxometry involve realizations of decoherence-free Bloch states for nonlocal symmetry 
species that did not exist before the measurement. The eigen states are fully deter- 
mined by three temperature-independent parameters which are effectively measured: 
the energy difference between tautomer sublattices; the double-well asymmetry for 
proton oscillators; the delocalization degree of the wavefunctions. The spontaneous 
^ ' decay of Bloch states accounts for relaxometry data. On the other hand, static 

H : 

, states realized by elastic scattering account for diffraction data. We conclude that 

both quantum and classical physics hold at every temperature. 
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I. INTRODUCTION 



The crystal of benzoic acid (BA, CgHsCOOH)^"— and that of potassium hydrogen car- 
bonate (KHCOs)^""— are reference cases for thermal interconversion of cyclic dimers, namely 
(C6H5COOH)2 or (HCO3 )2, linked through bistable 0-H- ■ -O hydrogen bonds. It has been 
observed for a long time that the coexistence of two configurations, say "L" for 01 — H ■ ■ ■ 02 
and "-R" for 01 ■ ■ ■ H — 02, leads to centrosymmetric configurations, LL and RR, exclusive 
of the less stable entities, LR or RL. The ratio RR : LL changes smoothly over a broad 
range of temperature without any structural transition. BA is traditionally regarded as 
a prototype for a class of chemical reactions involving proton transfer and tautomerism, 
that is of fundamental importance in many fields across physics, chemistry and biology. 
Experimental and theoretical studies are expected to enlighten the quantum behavior of 
hydrogen bonds in complex environments, with particular emphasis on proton tunneling 
and on the emergence of the classical regime at elevated temperatures. Published works 
can be categorized as follows: (i) quantum, for vibrational spectra;— '^^'^^'i^ (ii) semiclassi- 
cal, for coarse-grained measurements that cannot resolve individual quantum levels, what 
is essentially the case for solid-state NMR or quasi-elastic neutron scattering (QENS);->i^ 
(iii) theoretical, for computed Born-Oppenheimer (BO) potentials.—"— Not surprisingly, dif- 
ferent approaches lead to severe confiicts of interpretation showing that there is so far no 
satisfactory interpretational framework. 

As a matter of fact, experiments conducted with different techniques have been inter- 
preted in contradictory ways. For example, whereas the ratio RR : LL is commonly 
attributed to disorder,-"^ neutron diffraction and vibrational spectra provide evidences 
of quantum correlations excluding disorder.—*^ Further conflicts appear when the relax- 
ation rate measured with NMR/QENS is modelled by stochastic jumps of rigid proton 
pairs, induced by a thermal bath.-ii^ By contrast, the quantum regime on the timescale ~ 
10"^^ — lO^^^s is attested to by vibrational spectra, from cryogenic to room temperature,— 
and symmetry related selection rules for Bloch states extended throughout the crystal ex- 
clude local dynamics. In addition, the tiny proton-proton coupling terms effectively observed 
exclude rigid pairs and the discrete density-of-states of the crystal is markedly different from 
an incoherent thermal bath.— liiii^ The point at issue is whether these conflicts arise because 
the fundamental nature of the crystal is beyond the reach of measurements, or whether the 
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ad hoc interpretative frameworks are not appropriate to recognizing common features under 
the hght of different techniques. 

Theoretical studies of the crystal of BA with computational quantum-chemistry are rather 
scarce.—*^ The most detailed works deal with the isolated dimer of formic acid (FAD).— 
Within the framework of the BO separation, it is postulated that every experiment can 
be represented by classical particles moving along definite trajectories. This "local realism" 
leads to the notions of "reaction-path" , "promoter mode" and "transition-state" for "proton- 
transfer", which are strictly alien to quantum physics.— In fact, such theoretical studies 
are hardly conclusive because the dominant mechanism is largely model-dependent: either 
stepwise for molecular dynamic simulations,—'^ or concerted, of course, for ab initio cal- 
culations with rigid pairs.^° Whatever the mechanism, quantum chemistry emphasizes that 
proton transfer should be promoted by a shortening of the O- ■ -O bonds, what is in conflict 
with the adiabatic separation of protons and heavy atoms attested to by experiments.— "iii^^ 
It is, therefore, an open question whether or not computational chemistry can capture the 
essential features of tautomerism. 

Our purpose is to show that we can eliminate every conflict of interpretation if we adopt 
the radical viewpoint that the crystal is a macroscopic-scale quantum system, essentially 
decoherence-free, at every temperature below melting or decomposition, and on timescales 
ranging from vibrational spectroscopy to diffraction.— liiii^ii^ Consequently, the classical be- 
havior of nuclei and the intuitive notions forwarded by computational chemistry are aban- 
doned. This quantum view of a macroscopic-scale open system is clearly in conflict with 
the theory of decoherence stipulating an irreversible transition to the classical regime due to 
the exceedingly fast loss of quantum coherence leaking out into the environment, at a rate 
growing exponentially with the number of particles. On the contrary, we shall show that 
the sublattice of bonding protons is immune to irreversible decoherence and that this immu- 
nity increases with the size of the crystal. This quantum theory was first elaborated for the 
crystal of KHCO3 and confronted successfully with observations.— The key point is that 
bonding protons behave as fermions, what is, to the best of our knowledge, unprecedented 
in the context of hydrogen bonding. Our purpose is to extend this theoretical framework to 
the more complex crystal of BA and to re-examine the conflicts of interpretation engendered 
by classical physics under the light of the theory of quantum measurements. 

Our rationale is based upon fundamental laws of quantum mechanics and experimental 
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FIG. 1. Sketch of tautomerism (I) and proton interconversion (II). 

facts for BA. The vibrational spectra, like those of KHCO3, are consistent with adiabatic 
separation and asymmetric double- wells for the OH stretch modes (i/OH).— lii^i^ In addition, 
the spectra reveal a phenomenon that is not observed for KHCO3: the thermal population 
of a state at hvi = (54 ±6) cm~^ is correlated with the interchange of single and double CO- 
bonds (tautomerism), occurring without any change of the proton states.— This is consistent 
with markedly different timescales for tautomerism (Fig. [H I) and adiabatic interconversion 
of protons (Fig. [H II), hereunder referred to as "interconversion" for short. By analogy 
with the popular "phonon-induced tautomerism",™ it was conjectured that hui could be a 
phonon such that the ground state should correspond to the preferred tautomer at low 
temperatures and the excited states to a mixture of both tautomers sharing the same z/OH 
states.— 1^2. However, neutron diffraction studies by Wilson et al.- report a temperature law 
for the proton-site occupancy factor (H-SOF) that is at variance with both the proposed 
model and NMR/QENS data. For example, the estimated -R-SOF at 20 K is ~ 13% instead 
of ~ given by the other techniques. In order to establish our theory upon an unquestionable 
experimental basis, we have carried out new measurements with neutron diffraction of the 
crystal parameters at various temperatures. 

In Sec. mi we report the structural parameters measured between 6 and 293 K and we 
show that there is no statistical disorder for protons. The reported H-SOF's are significantly 
different from those of Wilson et al. and more accurate. They are consistent with other 
observations. In addition, the quasi-linear correlation with the mean CO-bond lengths is 
a new information. Van t'Hoff laws for protons and CO-bonds compare favourably with 
spectroscopic data. The quantum theory is presented in Sec. IIIII We insist that vibrational 
spectra are representative of double-well operators for nonlocal observables, which are fun- 
damentally different in nature from potential functions of classical coordinates. In IIII Al 
we set the Bloch wavefunctions for a lattice of bosons and we establish a new energy level 
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scheme consistent with the new diffraction data. In IIIIBI we show that antisymmetrization 
for a lattice of fermions yields a macroscopic-scale time-independent state with remarkable 
properties. In IIII CI we show that the crystal is immune to decoherence and we introduce 
the quantum-measurement scheme. Then, the theoretical temperature law for the i?-SOF 
is compared with measurements in IIIIDi In IIII El we show that NMR/QENS data are 
consistent with a coherent decay of measurement-induced transitory states. 

II. NEUTRON DIFFRACTION 

Data were collected with the four-circle diffractometer, 5C2, at the Orphee reactor (Lab- 
oratoire Leon Brillouin) .-Si Needle-shaped single-crystals (3 x 3 x 10 mm^) were obtained by 
slow cooling of an ethanol-hexane solution. A specimen wrapped in aluminum was loaded 
in a closed-cycle-refrigerator whose temperature was controlled to within ±1 K. The struc- 
tural parameters listed in Table [T] were obtained with CRYSTALS.— Positional and thermal 
parameters and the SOF's (p^, pn) were allowed to vary independently. 

A. Crystal structure 

Inspection of intensities for absent reflections confirms the P2i/c (C|;j) space group as- 
signment at every temperature. There is no evidence of disorder or domains. The unit 
cell parameters increase slightly with the temperature but there is no obvious discontinu- 
ity suggesting a phase transformation. There are four equivalent CeHsCOOH entities per 
unit cell. Centrosymmetric dimers are linked through hydrogen bonds whose lengths 2.61 
A < -Roo < 2.65 A are slightly longer than those of KHCO3 in the same temperature range 
(2.58 — 2.60 A). The mean distances between R- and L-sites are virtually temperature in- 
dependent: L- L = R- R ^ 2.24 A and L - ^ 0.70 A 0.60 A for KHCO3). The 
aromatic-ring is practically temperature- independent, the CC-bonds are equivalent and the 
angles equal 120°. 

Wilson et al. have pointed out that dimers in planar layers are connected to each other 
by contacts between oxygens and hydrogens of the benzene ring.- These H- ■ -O distances (a 
and b in Table [T]) increase slightly with temperature but there is no obvious correlation with 
SOF's or bond lengths. 
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B. Quantum order versus classical disorder 



The structural model proposed by Wilson et al., is based on conflicting premises: Frac- 
tional SOF's are supposed to be representative of a statistical distribution of mutually exclu- 
sive empty or occupied sites and, at the same time, they are constrained by pi(T) + pji{T) = 
1. This is simply incorrect because a random distribution of occupied and empty sites should 
transfer some Bragg-peak intensity to a background of coherent diffuse scattering leading 
to an apparent decrease of the number of protons effectively detected.— At elevated tem- 
peratures, this would give p^. ~ P_r ~ 0.25, instead of 0.50 measured. In fact, our own 
reflnements show that Pl{T) + Pr{T) converges to 1 without any constraint and fractional 
SOF's can be rationalized without recourse to disorder with the partial differential nuclear 
cross-section in the quantum regime:— 

d^ = ^EP^KWI(/inr)K)|k™)r5(;^a; + E,-E;). (1) 

' iJ 

The summation runs over initial states Ei, weight pi, and flnal states |/), Ef. The 
neutron wave vectors are k^j and k„j. The energy transfer is hu and dQ is an element of solid 
angle. The interaction operator, V"(r), comprises operators at every site. For diffraction, 
l^) = I/), Ei = Ef, hu = 0. Thanks to adiabatic separation, piiT) + Pr(T) = 1 is 
representative of distinct sublattices of LL or RR dimers, respectively. 

Alternatively, should the space periodicity be destroyed by disorder, the mean potential 
'^Pi{i\V{T)\i) could be Fourier transformed into the unit-cell structure factor and only the 

i 

mean position of protons would be relevant.— This model was graphically illustrated by 
Wilson et al.'^ who showed that protons located at L-sites at low temperature shift smoothly 
toward the center of the O- • -O bonds as the temperature increases, while the anisotropy 
of the thermal ellipsoid increases tremendously along O- ■ -O, as if the local potential were 
dramatically flattened. In contrast, the vibrational spectra do not show any signiflcant 
frequency shift. This model is untenable and so is statistical disorder. 

C. Temperature effects 

Fig. [2] shows substantial differences between Pr(T) in Table [T] (solid circles) and the 
i?-SOF's of Wilson et al. (X),- especially at low temperatures. Our estimate Pr{T) = 
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FIG. 2. (Color on line) Temperature laws for the occupancy factor of the less favored proton 
sites. • with error bars: experimental points (this work), x: experimental points after Ref. Q]- 
Blue dot dashed with bars: SOF due to tautomerism only. Solid plus grey zone: Tautomerism 
plus interconversion without the intermediate states (jl4p . Red dashed with bars: Tautomerism 
plus interconversion including the intermediate states (|12p . The grey zone and the bars represent 
dispersions of the theoretical curves due to uncertainties for the parameters. Curves A and B are 
analogous to (fH|) and (fT2|) . respectively, for a single tautomer state (see text). 

at 6 K is in accordance with a well defined minimum of the potential energy. There is no 
more conflict with other techniques. The best fit to our data is consistent with a van t'Hoff 
equation substantially different from that of Wilson et al.: 

In — = -0.07 ± 0.03 + (86.8 ± 1.7)T~^ (2) 
Pl 

The difference in enthalpy Aifp = (86.8 ±1.7) K compares with the NMR-based estimate: 
Ai7^Mij = (86.5 ±1.5) K.^ 

Fig. [3] suggests a linear correlation of the CO bond lengths and p^j, apart from the points 
at 100 K [pr ^ 0.30) for which the bonds are longer than expected by ^ 0.01 A, what is 
about twice the variance. Since at 100 K in Fig. |2]is not singular, there is no reason to 
question the data acquisition. Statistical errors could be actually greater than the variance 
given by CRYSTALS. 

Because we know from spectroscopy that the frequencies of the COO modes are temperature- 
independent,— the increase of the C=0 length and the decrease of the C-0 bond do not 
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FIG. 3. CO bond lengths as a function of the occupancy of the less favored proton sites. ■ and ♦ 
with error bars: experimental points (this work). □ and 0: after Ref. [3]. The straight lines are 
guides for the eyes. 

correlate with a continuous change from the carboxylic groupings COOH to the carboxylate 
COO^H+. The bond lengths are more likely representative of an unresolved mixture of 
tautomer sub lattices with constant CO-bond lengths. Then, with Rc=o = 1-233 A and 
Rc-oH = 1.320 A (Table [U 6 K), best fit exercises give: 



In ^^-""^ Z.,^^'^'^^^ = 0.01 ± 0.10 + (85.1 ± 6A)T~\ 
Rcioi\T) — Rc=o 

In ^^=^ --^^^f'^(^) = 0.18 ± 0.15 + (71.7 ± 9.5)T-\ 

Rci02\T) — Rc-OH 



(3) 



As these equations should represent the same equilibrium, we take the mean values ASco = 
0.1±0.2, AHco = (78±20) K, and the measured bond lengths are represented by a two-level 
system as: 

^cioilT") = Po{T)Rc=o + PiiT) Rc-OH, 

(4) 

Rci02{T) = po{T)Rc-oR + Pi{T)Rc=o; 
with po{T) = (1 + e-^^co/Tyi g^^^ ^^(2") = po{T)e-^"^o/^ . The solid lines in Fig. S 
and the dispersion zone (grey) representative of statistical errors compare reasonably with 
the experimental points. The mean deviation is less than 5 x 10"^ A for each branch. 
Alternatively, a harmonic energy level scheme for a phonon such that the excited states be 
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comprised of both tautomers in equal proportions,— ^iiii^ should give: 

i?cioi(T) = p'oiT)Rc=o + p'iiT)iRc=o + Rc-OR)/2; 
^ci02(T) = p',{T)Rc-ou+p[iT)iRc=o + Rc-ou)/2; 

with p[)(T) = (1 - e-^^co/T)-i g^j^^ = e-^^c:o/T_ rpj^g dashed curves in Fig. His 

less satisfactory. In fact, it is unlikely that a single phonon could drive at once the many 
degrees of freedom involved in tautomerism. In addition, a flip-flop of the carboxylic group 
is sterically forbidden. Presumably, the monoclinic crystal field separates the ground states 
of the otherwise degenerate tautomer sublattices and AHqo is not representative of an 
asymmetric double-well along a definite reaction path. 




1 1 1 1 1 1 1 1 1 1 — 
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FIG. 4. CO bond lengths as a function of the temperature. ■ and ♦ with error bars: experimental 
points from Table HI Solid curves were calculated with ^ and Rc=o = 1-230(5) A, Rq-ou = 
1.318(5) A, AE = (78 it 20) K~^. The grey zones represent the statistical dispersion. The dashed 



curves ([5]) hold for the coexistence of both tautomers at elevated temperature 



Finahy, the various estimates AHqo, ^Hp, AHnmr, and hui = (81 ± 9) K,—^ yield the 
mean value of AH = (78 ± 20) K. In Fig. [21 the dot-dashed curve pn = Pi(T) is slightly 
underestimated above ~ 150 K, what suggests that adiabatic interconversion (Fig. [H II) 
is not negligible. The purpose of the quantum theory presented below is to calculate this 
contribution on a physically sound ground. 



III. QUANTUM THEORY 

The stretching modes play the leading role for adiabatic interconversion. In classical 
physics, a double- well between L1L2 and R1R2 is permitted for the symmetric coordinate Xg 
(Raman- active), while it is excluded for the antisymmetric Xu (infrared-active). The effective 
mass fig = 2 amu is eventually increased by a promoter mode involving heavy nuclei.— 
However, knowing that classical trajectories are alien to quantum mechanics, vibrational 
spectra effectively reveal the eigenstates of observables, which are not classical coordinates. 
For a periodic lattice of indistinguishable centrosymmetric proton sites, the observables are 
symmetry species at the center of the Brillouin-zone (CBZ), say X^j^CjC^) with C, = "tt" or 
"(7", for the stretching, in-plane bending, out-of-plane bending modes, respectively, and it 
is an experimental fact that both infrared and Raman spectra are consistent with double- 
well operators V^{xi)y and single- well operators V^{ri^), V^{Q): ^oi proton oscillators whose 
effective masses are /i^ = /i^ = 1 amu (see Table HB . The INS spectrum is also consistent 
with this interpretation and reveal negligible dispersion for the OH modes. The eigenstates 
are, therefore, entirely determined. 



A. The OH states 

For the ease of notation, an ideal disorder-free crystal of BA can be represented by A/" 
reduced unit-cells indexed j in 3-dimension, each of them comprising only one planar dimer. 
JV is on the order of Avogadro's constant. With respect to the center of symmetry, the mean 
positional parameters of the bonding protons are {xoj,yoj,0} for Lij, {— xqj, —yoj,0} for 
L2,j, {—Xoj,yoj,0} for Rij, {xqj, —yoj,0} for R2J. Small displacements are represented 
by local (classical) variables sij,S2j (s = x,y,z), and the mass- and distance-conserving 
normal coordinates are: 

1 , . 1 , 

^Lu,j = ^(si,i + S2,i -2soj); SL„j = — =(sij-S2j); 

_ Y _ r (6) 

^Ru,j - y^l-^ij + -52^ + 2so,j); sr^j - -^(sij - S2,j); 

In classical physics, the effective oscillator mass of a normal mode is arbitrary. In quan- 
tum physics, normal coordinates are observables and their effective masses are definite and 
measurable.— 
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For each symmetry species, there are M Bloch states indexed r (1 < r < M) whose 
wavefunctions are hnear combinations of those of the unit-cells. For the lower states of the 
double- well V^(xg), they can be written as: 

M 



VAT i=l 
1 ^ 

'0i^r(kir,t, V^ir) = -7=y^[-sin 9^ ipLU + cos9^ V^/j^j] expi[kir.Hj + wi^^t + V^i^r] 



(7) 



\^nr is a wave vector; Hj is a vector of the sublattice; (fn^r is an arbitrary phase.^- As there 
is no permanent electric dipole for centrosymmetric dimers, dispersion is negligible and the 
states are degenerate with respect to k„r- Each wavefunction spans every proton site and 
the occupation factor at a given site includes contributions from every state. ipL£,,j and ipR^^j 
are harmonic eigenfunctions in 3-D for second-order expansions of the potential operator 
y^(Xg) + y^iVi) + ^(Cc) around each minimum. Even at 300 K, only the ground states of 
the bending modes are significantly populated and the corresponding quantum numbers are 
omitted for the ease of notation. Because of the double-well, the symmetry species Xi are not 
linear combinations of normal coordinates (|6]) and, conversely, ([7]) contains no information 
on the local coordinates Xij,X2^j. 

The eigenstates of the z/OH modes given in Table [III are consistent with: hv^i^ = 
hi^oig = huoi, = Og = 9, tan26' = I'tJ {^oi — ^t), where hut = (6 ± 1) cm~"^ is the 
tunnel splitting calculated numerically for the symmetric double- wells. Then, cos ^ ^ 1 and 
sin^ = 6 = (1.8 ± 0.3)10^^.^^ The lack of u — g splitting for huoi, as opposed to higher 
transitions, is a consequence of the marked localization of the wavefunctions in either well. 
Creation(annihilation) of a quantum huQi can be, therefore, conceived of as the transfer of a 
particle-like entity from one minimum L^{R^) to the other one R^{L^), without any energy 
transfer to the internal degree of freedom. 

In classical physics, interconversion is thought of as an equilibrium between dimers L1L2 
and R1R2, knowing that L1R2 and R1L2 are energetically excluded. Each dimer entity 
is supposed to be at any given time in a definite configuration, whether it is observed or 
not (classical realism). In quantum physics, symmetry species yield a much richer 4-level 
scheme. For the preferred tautomer sublattice at low temperatures (say To, labeled with 
the quantum number / = 0), the four-well operator VoiXujyXgj) = ^OuiXuj) + ^OgiXgj) 
is shown in Fig. |5l The central barrier of ~ lO'* cm~^ is twice the barrier in 1-D. It is 
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FIG. 5. Perspective view of the potential operator for the hydrogen bonding protons in the crystal 
of benzoic acid. Xu and Xg the symmetry species of the OH stretching mode unveiled by 
vibrational spectra. 

surrounded by four deep minima: L^Lg set to zero-energy; RuLg and L^Rg at ^ huQi] 
RuRg, at ~ 2huQi. Interconversion between LuLg and RuRg occurs through RuLg and LuRg 
at a very low energy cost compared to the classical mechanism, via R1L2, L1R2. For the less 
favored tautomer sublattice (Ti, / = 1), the C— O and C=0 bonds are interchanged. The 
operator Vi{Xu,j, Xg,j) = ViuiXuj) + ^igiXgj) is obtained by changing the signs of the linear 
terms, while the other terms are unaffected. Consequently: RuRg is more stable that LuLg] 
([7]) transforms as 9 — )■ n/2 — 9; huoi, huos/., are unchanged. 

The scheme presented in Fig. [6] is consistent with the diffraction data reported in Sec. [Ill 
(Fig. H]). Compared to that previously published,— the main difference is that the mixture 
of To and Ti for every Z > 1 is now replaced by a two-level scheme: / = for To and I = 1 
for Ti. The wavefunctions are: 

1 ^ 

J'lnungi^nurigry'ty '^InuTigr) = j 6Xp j [k„^„^r -Hj + iOln^Ugt + '^In^rigr]'-, (8) 

i 

where ^^zn^n^j = i^inu,j^ing,j-— The probability density l^'ooojP, or l^oiijP, is localized to 
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FIG. 6. Schematic view of the wave functions of the unit cell j for the hydrogen bonding protons 
in the crystal of benzoic acid. Xuj and Xgj are the symmetry species of the unit cell unveiled 
by vibrational spectra. For the sake of clarity, the weak component of the wave function in one- 
dimension is multiplied by a factor of 5. The quantum numbers are l,nu,ng, respectively. 



order around L^Lg, or RuRg, respectively. The mean geometry is L1L2 or -Rii?2- By 
contrast, |^I/oio,iP a^nd I^I/ooijP have no counterpart in classical physics. They represent 
counterintuitive configurations with equal occupation factors at each of the four sites. 
This is markedly different from the transition state of computational chemistry, with 
two protons at the center of the O- • -O bonds. For / = 1, energy levels are shifted upwards 
by hui and the wave functions are sketched in Fig. [6]-rhs. 



13 



B. The crystal states 



Bloch states obey the Bose statistics law. If the fermion character of protons can be 
ignored, the sublattice can be represented by a superposition of degenerate states: 

)); (9) 

l,nu,ng r 

«wo = + +Poi)'^; alio = «m = «woPoi; afn = aLPoi; Pi = exp{-hui/kT); 

Poi = exp(— /iz/Qi/ZcT). Then, destructive interferences cancel the phase factor, in space and 
time, so the probability density is zero everywhere, at any time, the Bloch states do not 
exist and the sublattice of protons is indefinite. 

Alternatively, when the fermion character of bare protons emerges from both the ionic 
character of the hydrogen bond and the adiabatic separation, Bloch states are forbidden 
by the exclusion principle. For a lattice of noninteracting bare protons (no spin-spin or 
electric-charge interaction, no exchange or Fermi energy) protons behave as sharply localized 
particles whose wavefunctions, say $(r— roj) at a given site in a given unit cell, do not overlap: 
J <|)(r — roj/)$(r — roj)(ir = 6jj'. This yields zero-width eigenf unctions, say Cy, C2j, T^ij, '^2j, 
for which mean-square amplitudes equal zero and tunneling is excluded. The sublattice of 
fermions is rigorously static and periodic (distortion-free). The centrosymmetric dimers 
CijC2j and 7lij7l2j correspond to eigenstates \00)j and respectively, such that Em — 

Em = (— l)'2/ii/oi- The states \7lijC2j)j and \CijlZ2j)j are excluded. There is no further 
eigenstate and the potential operator is irrelevant. In addition, the antisymmetry postulate 
for centrosymmetric pairs requires spin state- vectors with singlet-like antisymmetry, \S)j, 
or triplet-like symmetry, |T)j,— such as: 

= ^ ^S'2ij[(aooo + ttiii)|00)j (aoii + aioo)|ll)j]|'5)j 

1 . (10) 
+ -j=2^A2ij[{aQ0Q + q;iii)|00)j + (aoii + aioo)|ll)j]|7')j; 

^21] = (1 + -P2ii)/2 is the symmetrizer; A2ij = (1 — p2ij)/2 is the antisymmetrizer; p2\j 
is the permuter;^^ af^Q = p\{l +pi)~^{l +Poi)~^ ^ni — '^roPoi- The thermal statistics are, 
therefore, different for fermions and bosons. However, energy transfer to the static lattice 
and interconversion between CijC2j and TZijTZ2j are forbidden and it is necessary to suppose 
a thermally induced mechanism to account for thermal equilibration. 
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C. Quantum measurements and decoherence 

In the first step of a quantum measurement, the system to be measured is entangled 
with another quantum object so as to reahze, or prepare, a state that will be effectively 
measured. At the microscopic level of a single object, the back-action of the measurement 
yields unavoidable changes of the system due to the collapse of the wavefunction. In classical 
physics, by contrast, a measurement can be carried out with unlimited precision without 
affecting the system being measured (noninvasive measurement). From the viewpoint of 
measurements, a crystal of BA is both quantum and classical. First, an incoming wave 
can realize a pure state (r). Second, the pre-existing state is practically unaffected by the 
realization of a state amongst a quasi-infinite number of equally probable realizations. The 
back-action can be ignored. In addition, the remote detection of the outgoing wave has no 
back-action to the measured system. Third, the probability for an incident wave to interact 
with a previously realized state is insignificant. The evolution of the prepared state is deter- 
mined by the time-dependent Schrodinger equation, irrespective of further measurements. 
Fourth, both the pre-existing and the reahzed states depend on the temperature that is a 
macroscopic-scale classical variable of the crystal. 

From the viewpoint of decoherence, the crystal is an open system in a gaseous atmosphere. 
Under normal conditions, the density of the surroundings is about 1000-fold less than that 
of the crystal, so the probability for an incoming thermal photon to realize a proton state is 
significant. However, decoherence is hampered because quantum correlations are constrained 
by the crystal structure. For bosons (Q, creation/anihilation of a Bloch state does not induce 
decoherence. For fermions ffTOl) . there is noway to transfer any energy to the static state, so 
as to reach thermal equilibrium, other than via realization of a driven Bloch state violating 
the quantum laws for fermions. We show below that the lifetime of such states is short 
enough to avoid irreversible decoherence of the bulk. 

In order to confront the theory with observations, we examine which experiments yield 
different outcomes for fermions or bosons. In the first instance, we can discard vibrational 
spectroscopy because the same Bloch states are realized via resonant energy transfer in 
both cases. By contrast, neutron diffraction and relaxometry provide experimental support 
to fermions. 
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D. Neutron diffraction 



Elastic momentum transfer (Q = k„j — k„/) realizes a Bloch state whose wave vector is 
Q. For uncorrelated scattering events, the state spans the M unit cells, irrespective of the 
coherence length of the beam. In the boson case, the realized state is 

\Bl)r= ^ a/n„nJ^«n„ng(Q,t))- (H) 

Diffraction gives Bragg-peaks whose intensities are proportional to the nuclear coherent 
cross-section for protons (cTc ~ 1.8 b) and to the Debye- Waller factor, an exponentially 
decreasing function of Q due to distortions. The occupancy factors are L-SOF = a^^Q + 
(floio + «ooi + flno + «?oi)/2 + flni, ^-SOF = af^Q + {aj^^ + aj^^ + ag^o + alo^)/2 + ag,,, and 
L-SOF -|- i?-SOF = 1 at every temperature. The temperature law for i?-SOF 

"^^'^^ = (1+P,)(1+P0,)' 

gives the dashed curve in Fig. [2] that is overestimated compared to measurements. 

For a lattice of fermions, the static Bloch state effectively realized depends on Q. For 
Q = L* 7^ H* (L* is the reciprocal crystal lattice), Q-induced distortions destroy the spin- 
symmetry: 

\Ht)r = --^^[(aooo + aiii)|00)j-h (aoii + aioo)|ll)j]exp(zQ.Lj). ^i^-^ 

The diffraction pattern is basically the same as that anticipated for bosons. However, the 
occupancy factors are different: L-SOF = agoo + Q^iii! -R-SOF = af^Q + alu] L-SOF -|- 
i?-SOF = 1. The temperature law for i?-SOF, 

PRFiT) = ,y (14) 

(1 +Pi)(l +Poi) 

gives the solid curve and the grey dispersion zone (Fig. [2]) in better agreement with experi- 
ments than Prb(T).— 

For the sake of comparison with KHCO3, the curves A and B in Fig. [2] represent (IT^ and 
f[T^ . respectively, for a single tautomer {hui ^ kT). These curves are markedly different and, 
for KHCO3, measurements are unambiguously in accordance with fermions fll4p . whereas 
bosons (fT2|) are excluded.— For BA, the difference is not so dramatic, because of the 
tautomer statistics, but the same trend in favour of fermions is observed. 
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It should be mentioned that for Q = H* the Bloch state \H2)r = |F)/VAr retains the 
essential features of the pre-existing state: static periodicity and spin-symmetry. Then, the 
cross-section for coherent scattering is the total cross-section including the huge incoherent 
scattering cross-section (cTj ~ 80 b): a = cXi+cXc ~ 82.0 b.^^ The enhanced diffraction pattern 
H* superimposed to the regular pattern L* is an evidence of the spin-symmetry featuring 
fermions. In addition, because the Debye- Waller factor equals 1 at every temperature, 
the contrast is best at large Q-values and at elevated temperatures. This was effectively 
observed for KHCO3 and, as expected, not for KDCO3. So far, there is no such experiment 
ever reported for BA. 

E. Relaxometry 

In classical physics, relaxometry techniques are thought of as noninvasive measurements 
of predetermined fluctuations.-'^'^ QENS is specific to protons, for incoherent scattering by 
other nuclei is negligible. This is not so stark for NMR because the spin-lattice relaxation rate 
(T^^) depends on many interactions. Nevertheless, both techniques give for BA comparable 
relaxation rates, which are essentially due to the bonding protons.^ It is commonly supposed 
that uncorrelated jumps of proton pairs between L1L2 and -R1-R2 can be represented by a 
rate equation for the probability PL{t) that, given a LiL2-dimer at t = 0, the same dimer is 
R1R2 at time t, and conversely for PR{t). The stochastic process occurring over a relatively 
long timescale gives rise to a very weak Lorentzian profile in u, centered at a; = 0, whose 
width is proportional to the inverse relaxation time ^--^J^M-^ This scheme is clearly in 
conflict with the rule stipulating that the interaction of a wave and a crystal is nonlocal. 

In quantum physics, an incoming photon or neutron can realize a Bloch state of the 
stretching mode, whose wave vector is the momentum transfer along the hydrogen bond 
{Qx ~ for NMR), while the bending modes can be ignored. For a lattice of bosons, the 
realized state analogous to (fTTI) 

\B2)r = S 0.lnung\J^lnung{Qx,t)) (15) 

is stationary and does not contribute to any relaxation mechanism, thanks to adiabatic 
separation. The rate effectively observed is a prima facie evidence that protons do not 
behave like bosons. For fermions, off-resonance energy transfer [huj ^ hvi, hv^i . . .) realizes 
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a time-dependent state analogous to (fT5|) and sets the starting time t = for subsequent 
decay. In order to show that this decay can account for the observed rate, we consider, in 
the first instance, a simple superposition, say aooo|000) + aoio|010) (from now on, we ignore 
r and {LnuUg) = \ J^inung{Qxi t)))- The probability density in the stationary regime is to order 



|aooo|000) + aoio|010)|2 

(16) 

Oooo^^oo^V'loOg + «oio^«oi„V'|oOg + aoooaoio£(V'L - '^lJ^l, cos27rz/oit. 
Harmonic oscillations at frequency z/qi of si'ipR^ ~ from e%p']^/2 to in one well 

and from —£il)\^/2 to eil)\^/2 in the other well, during a half period (2z/oi)~^, determine the 
inverse residence time: 2£Vqi = 2//^, in proton per second units (Hs^^). (Similar terms are 
given in Table UTTl for every binary combination.) Then, we define the relaxation rate as the 

1 1 /2 

critical decay of the normalized amplitude: Tq^ = 2£:z/oiaoio/'^ooo = '^^^oiPoi , and likewise 
for TQg^. The fractional exponent is distinctive of a coherent process. 

Next, we consider the probability density | X]n„ng '^"u'^g 1^^"%) P ^'^^ ^ single tautomer 
{Pi = 0)- Because cannot distinguish u and g species, we obtain two rates (to order e^): 
|aoo|00)+aoi|01)+aio|10)p gives Tq"^ = Tq^^+Tq^^ = 4£z/oipJf and |aoi|01)+aio|10)+aii|ll)|2 
gives 

-1 -1,-1 A CtOlOll . 3/2 

"-00 

In fact, because the four states are correlated, the coherent decay limited by the longer res- 
idence time Ti > To gives the overall rate 2t^^ = 8ep'^{^. This accords with the temperature 
law effectively measured for KIICO3.— 

For BA, NMR and QENS are both consistent with a markedly different empirical law:- 

T-l = (1.72 ± 0.02)108 coth[(43 ± l)/kT] ^^^^ 
+ 10^0 exp(-120/A;T) + 6.3 x lO^^ exp(-400/A;T). 

The first term was ascribed to the transition rate of the thermal bath supposed to rock the 
double-well asymmetry of (86 ± 2) K for LL and RR. The other terms were attributed to 
thermal jumps in the classical regime. Alternatively, Fig. [7] shows that (fTSl) (solid curve) 
is consistent with the theoretical law for coherent decay of a Bloch state (dashed curve and 
grey zone): 

Tq' = 8eW coth{hui/2kT) + Ae'^mpH' [p/"'^' + p'/'] + ^empli^- (19) 
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FIG. 7. Comparison of the experimental temperature law of BA (solid line) (jlSp with the quantum 
theory (dot-dashed) ([IS]): uqi = 5.16 x lO^^g-i^ ^ ^ g.OlS, r^^ = (2.3 ± 1)10^ coth[(39 ± 10)/r] + 
(6.77 ± 1.15)109{exp[-(89 ± 13)/r] + exp[-(167 ± 13)/r]} + (7.2 ± 1.2)10^1 exp[-(385 ± 9)/r], 
in proton per second units. The grey zone features the dispersion of the theoretical curve due to 
uncertainties. The white curve was obtained with e = 0.016. 

This function is different from that derived from the previous energy level scheme.— It 
contains more terms than f|T8|) but it is more economic, since six empirical parameters are 
replaced by three scaling factors {ui, z/oi,e) bridging relaxometry and spectroscopy. 

The first theoretical term is formally identical to that in (fT8|) but it is fundamentally 
different in nature. The rationale is that thermal transitions at ±/iz// occurring on a much 
shorter timescale than the interchange of single and double CO-bonds do not toggle the 
double-well asymmetry. The cross-section for one-quantum annihilation is the Bose factor 
n(z/;) = [explhui / kT) — the cross-section for one-quantum creation is [l-|-n(z/i)] and the 
population ratio is n(z/;)/[l + nlui)] = pi?^ Annihilation yields n{vi) tautomer lattices Ti in 
the lower energy level (/ = 0) and creation yields [1 + nijJi)] To in the upper level (/ = 1). 
The total population of Tq and Ti occupying the same level is [1 -|- 2?2(z/;)] = coth(/iz/;/2A;T), 
that is the mean energy in hvi units of the two-level system. Then, an incoming wave can 
realize 



H^) = coth{hui/2kT)\B^). 



(20) 
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The superposition of proton eigenstates for To and %i occupying the same level, 

coth{hui/2kT) ^ a„„„J|On„ng) + |ln„ng)], (21) 

riurig 

gives rise to a tunnel splitting proportional to the overlap integral, /iz/± = hh't{Onung\lnung) = 
Ae'^i'oi, and the decay rate is 2hi'± coi\i{hvi/2kT). We are not aware of any published exam- 
ple of a tunnel splitting promoted by thermal fluctuations, or, at very low temperatures, by 
quantum fluctuations of the zero-point energy. It is rather counterintuitive that a coherent 
superposition of proton states could emerge from an incoherent mixture of tautomer states. 
This is possible for fermions because the pre-existing state is insulated from environmental 
fluctuations and the coherent superposition of the four-well operators for Tq and %i is real- 
ized at the time of the measurement for a seemingly static mixture of tautomers in the same 
state. For bosons, the same fluctuations would yield a mixture of the pre-existing four-well 
operators and no tunneling. 

The second term in square brackets corresponds to combinations such as aooo|000) + 
aiio|110) + aioi|101) and aioo|000) + aoio|010) +aooi|001), which should not occur for bosons. 
The rate is independent of coi\i{hvi / 2kT) because it does not matter which is occupying 
the lower/upper level. The longest residence time is (2e^i/oi)~^ (see Table UTTl) and the rate is 
proportional to ano/aooo = aioi/aooo = vixVi'^ and aoio/aioo = ctooi/aioo = vixVl^^'^ ■ The 
composite nature of this term is overlooked in (fT8|) but Fig. [7] shows that the theoretical 
curve is close to observations in the intermediate temperature range (50 — 150 K). Finally, 
the last term of (fT9|) is identical to (fT7|) . 

The coherent decay of the real part of the wavefunction fl20|) accounts for the observed rate 
at every temperature. (The white curve in Fig. [7]shows that the agreement is improved with 
e = 1.6 X 10~^.) This deterministic evolution is not related to "proton transfer" in any way. 
The scaling factors p]^'^ ,pl{^ ,e, are effectively measured and numerics compare favorably 
with spectroscopic data. If this agreement is not fortuitous, environmental decoherence, 
transition rate induced by a thermal bath and semiclassical jumps should be excluded. In 
addition, the imaginary part of the wavefunction is fully determined by the measurement 
setting Qx and t = 0. We can thus obtain on a purely experimental basis the maximal 
knowledge that is theoretically available about the state. The lifetime of the Bloch state 
(~ 10^'' — lO^^^s) is much longer than the period = 2 x lO^^^s of the z/OH mode and 
quasi-stationary quantum interferences occur on a macroscopic scale at every temperature. 

20 



Interesting enough, the sphtting hi>± 8 x 10 ^ cm ^ estimated for the crystal of BA 
compares favorably with those reported for isolated dimers in supersonic jet expansions 
for which hui = 0: ^ 3 x lO^^ - 16 x 10"^ cm^^ for (DC00H)2; ^ 12 x lO^^ cm-^ for 
(HCCOH)2;-2S ^ 8x lO^^ cm-^ for formic-acetic acid (HCOOH...HOOCCH3);37 ^ IQ-^ cm-^ 
for formic-propiolic acid (HCOOH...HOOCCCH).— In these cases, the splitting should exist 
whether it is observed or not. The four- well operator of BA is apparently a feature of 
carboxylic-acid dimers, centrosymmetric or not, irrespective of the physical state. It is also 
rather counterintuitive that isolated dimers at very low temperature compare with BA in 
the high temperature regime: kT ^ hui. 

IV. CONCLUSION 

The structure of the crystal of benzoic acid is consistent with a bistable system of unre- 
solved tautomer sublattices sharing the same Hydrogen-bonding proton-sites and there is no 
evidence of disorder for these protons. Diffraction, vibrational spectroscopy and relaxometry 
techniques, provide evidences of wave-like protons excluding definite particles. The notions 
of reaction path across a potential function, transition state, stepwise or concerted proton 
transfer, are irrelevant. 

We propose a quantum theory for the bonding protons regarded as fermions. The exclu- 
sion principle and the antisymmetry postulate yield a superposition of static lattice-states 
with spin-symmetry. According to the theory of quantum measurements, an incoming wave 
is supposed to realize a pure state without any significant back-action to the lattice state. 
For vibrational spectroscopy, resonant energy transfer to the pre-existing static lattice re- 
alizes the time-dependent eigenstates of a four-well operator in the abstract space of the 
symmetry species. For NMR and QENS, off- resonance energy transfer realizes a transitory 
state and the coherent decay rate of the real part of the wavefunction is measured. For 
diffraction, elastic momentum transfer realizes a static proton state in accordance with the 
site occupancies effectively measured. The theory accounts for every measurement with 
only three temperature-independent factors. Observations supporting a pre-existing lattice 
of fermion are the proton-site occupation factors revealed by diffraction, on the one hand, 
and, on the other, the spontaneous decay of Bloch states and the tunnel splitting promoted 
by environmental fluctuations. The alternative possibility that protons could behave like 
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bosons when they are not observed can be ruled out. 

Both quantum and classical physics hold for the nuclei of an open crystal of BA at 
every temperature. Macroscopic-scale quantum correlations rooted in the periodic structure 
cannot be destroyed and there is no dividing line between the quantum and the classical 
regime. 
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TABLE I. Neutron single crystal diffraction data and structure refinement. A = 0.8305 A. Space group P2i/c. a{l) limit: 3.00. Refinement 
on F. The occupation numbers for the proton sites are pL and pR. The variance for the last digit is given in parentheses. 



Crystal data 


6 K 


25 K 


50 K 


75 K 


100 K 


125 K 


150 K 


175 K 


200 K 


250 K 


293 K 


a(A) 


5.401(1) 


5.401(1) 


5.415(1) 


5.415(1) 


5.500(1) 


5.427(1) 


5.456(1) 


5.456(1) 


5.500(1) 


5.500(1) 


5.500(1) 


6(A) 


5.004(1) 


5.004(1) 


5.016(1) 


5.016(1) 


5.100(1) 


5.112(1) 


5.074(1) 


5.074(1) 


5.100(1) 


5.100(1) 


5.100(1) 


c(A) 


21.879(1) 21.879(1) 


21.826(1) 


21.826(1) 


22.020(1) 21.810(1) 21.875(1) 21.875(1) 22.020(1) 


22.020(1) 22.020(1) 


P 


98.47(1)° 


98.47(1)° 


98.44(1)° 


98.44(1)° 


97.90(1)° 


98.26(1)° 


98.12(1)° 


98.12(1)° 


97.90(1)° 


97.90(1)° 


97.90(1)° 


Volume (A^) 


584.8(2) 


584.8(2) 


586.4(2) 


586.4(2) 


611.8(2) 


598.8(2) 


599.4(2) 


599.4(2) 


611.8(2) 


611.8(2) 


611.8(2) 


Reflections measured 


2950 


3594 


2797 


3364 


3067 


3092 


2778 


3607 


3502 


3828 


3401 


Independent reflections 


2589 


2637 


1946 


2645 


2724 


2283 


2297 


2300 


2714 


2748 


2630 


Reflections used 


1649 


1679 


1382 


1472 


1363 


1122 


1018 


951 


1031 


755 


511 


R-factor 


0.093 


0.079 


0.058 


0.090 


0.051 


0.048 


0.047 


0.041 


0.069 


0.049 


0.50 


Weighted R-factor 


0.114 


0.097 


0.068 


0.051 


0.053 


0.048 


0.049 


0.037 


0.045 


0.032 


0.030 


Number of parameters 


132 


137 


137 


137 


137 


137 


137 


137 


137 


137 


146 


Goodness of fit 


1.148 


1.117 


1.113 


1.026 


1.085 


1.124 


1.134 


1.101 


1.033 


1.113 


1.101 


Occupation pL 


1.00 


0.966(13) 


0.852(13) 


0.761(13) 


0.69(1) 


0.65(1) 


0.61(1) 


0.59(1) 


0.583(12) 


0.57(2) 


0.50(1) 


Occupation pR 


0.00 


0.034(13) 


0.148(13) 


0.239(13) 


0.31(1) 


0.34(1) 


0.38(1) 


0.41(1) 


0.417(12) 


0.43(2) 


0.50(1) 


diR (A) 






0.68(1) 


0.71(1) 


0.69(1) 


0.69(1) 


0.67(1) 


0.72(1) 


0.70(1) 


0.69(2) 


0.68(1) 


Roo (A) 


2.623(5) 


2.621(5) 


2.619(5) 


2.614(5) 


2.646 (5) 


2.628(5) 


2.621(5) 


2.610 (5) 


2.624(5) 


2.623(5) 


2.629(5) 


Bond length ClOl (A) 


1.233(3) 


1.236(3) 


1.244(3) 


1.251(3) 


1.268(2) 


1.261(3) 


1.269(3) 


1.263(2) 


1.271(3) 


1.264(3) 


1.260(6) 


Bond length CI 02 (A) 


1.320(3) 


1.313(3) 


1.305(3) 


1.293(3) 


1.302(1) 


1.286(3) 


1.279(3) 


1.280(3) 


1.279(3) 


1.276(4) 


1.266(6) 


contact a (A) 


2.546(3) 


2.545(3) 


2.531(3) 


2.530(3) 


2.549(4) 


2.522(3) 


2.536(3) 


2.541(4) 


2.565(3) 


2.563(5) 


2.610(6) 


contact h (A) 


2.462(3) 


2.466(3) 


2.480(3) 


2.490(3) 


2.536(4) 


2.517(3) 


2.528(3) 


2.534(4) 


2.548(3) 


2.553(8) 


2.553(6) 



TABLE II. Double-well operators, distances between minima (A^^), barrier heights (H^), and 
energy levels, for the crystal of the ring deuterated Benzoic acid, xs, ^re the symmetry species 

with inelastic neutron scattering after Ref. [9|]. 





V^ix^) (cm-i) 


Axo (A) (cm-i) 


huQ^ (cm ^) 


(cm ^) 


/ii/^3 (cm 1) 


Infrared ^'^ 


265x« + 0.2860 x lO'^Xu 


0.70 5005 


(172 + 4)* 


2570 + 20 


2840 + 20 




+171480 exp(-2.17xS) 










Ramar>ii 


270xg + 0.2829 x lO^Xg 


0.69 5006 


171 + 4 


2602 + 4 


2853 + 4 




+171120 exp(-2.15x3) 











TABLE III. Leading terms of the inverse residence time for binary superpositions of the dimer- 
oscillator states sketched in Fig. [6l 

|000) |010) |001) |011) |100) |110) |101) |111) 

|000) — 2ez^oi 2ei/oi 4:e'^'^oi Se^J^oi 2e^z^oi '^e'^'^oi 4e^foi 

|010) 2ei/oi — — 2euoi 2e^vn 8e^z^oi Se^^^oi ^e^t^oi 

|001) 2evQi — — 2ez^oi 2e^z^oi Se^'z^oi Se^i^oi 2e^foi 

|011) 4e^foi 2ez.'oi 2ei^oi — 4e^z^oi 2e^z^oi 2e^z^oi Se^^'oi 
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